Observation of spatial quantum correlations induced by multiple scattering of 

non-classical light 
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We present the experimental realization of spatial quantum correlations of photons that are in- 
duced by multiple scattering of squeezed light. The quantum correlation relates photons propagat- 
ing along two different light trajectories through the random medium and is infinite in range. Both 
positive and negative spatial quantum correlations are observed when varying the quantum state 
incident to the multiple scattering medium, and the magnitude of the correlations is controlled by 
the number of photons. The experimental results are in excellent agreement with recent theoretical 
proposals by implementing the full quantum model of multiple scattering. 
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Multiple scattering of light waves propagating through 
a disordered medium can be described as a random walk 
of the direction of propagation whereby the medium 
becomes opaque ^2] . Such random scattering determines 
the performance of nanophotonic devices based on, e.g., 
photonic crystals [3] or surface plasmon polaritons 4j, 
and can be applied for enhancing the information capac- 
ity for optical communication [5]. In multiple scatter- 
ing, the different possible light trajectories through the 
random medium will interfere and a complex spatial in- 
tensity distribution of light is generated, cf. Fig. [ijb). 
Usually interference effects tend to vanish after averag- 
ing over all possible realizations of the disorder and the 
light transport is described by diffusion theory. Under 
conditions where multiple scattering is strong, light dif- 
fusion is modified leading to spatial correlations such as 
mesoscopic fluctuations [1, 0, S] ■ These correlations im- 
ply that the light intensity observed at one position after 
the medium depends on the intensity at a different posi- 
tion even after averaging over all ensembles of disorder. 

The significance of the quantum nature of light in a 
multiple scattering context has only been addressed re- 
cently [l^, 11, l3, 13, Pioneering work has pro- 
vided the theoretical framework for a quantum optical 
description of multiple scattering that was used to 
predict the existence of spatial quantum correlations in- 
duced by multiple scattering of non-classical light 
These quantum correlations exist even in the diffusive 
regime of multiple scattering and add to the traditional 
mesoscopic correlations that are of classical origin Q. 
Here we present the experimental demonstration of spa- 
tial quantum correlations induced by multiple scattering 
of squeezed light. 

Due to the Heisenberg uncertainty principle intrinsic 
quantum fluctuations are always present. Although al- 
ways of inherent quantum origin, photon fluctuations 
are classified as either classical or non-classical if simi- 



lar fluctuations can be induced by classical light sources 
or not. The classical boundary corresponds to Poisso- 
nian photon statistics where the variance of the photon 
number fluctuations equals the mean number of photons. 
A purely quantum regime exists where photon fluctua- 
tions are reduced beyond the classical limit leading to 
sub-Poissonian photon statistics. 



Non-classical 
light source 




Spatial 
quantum 
correlations 




I Max 



I 



IVlin 



FIG. 1: (color online), (a) Illustration of different light tra- 
jectories in a multiple scattering medium that are correlated 
depending on the quantum state of the light incident to the 
medium, (b) Measured spatial intensity distribution of light 
transmitted through a multiple scattering medium displaying 
a volume speckle pattern. 

The number of photons exiting a multiple scattering 
medium in a specific direction can be correlated or anti- 
correlated with the number of photons in another direc- 
tion, and this spatial correlation depends on the quantum 
state of light illuminating the medium, cf. Fig[lja). The 
degree of correlation is quantified by the spatial quantum 
correlation function lisll 
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where a labels the direction of the incident light and h 
and h' the two different output directions. The operator 
flab represents the number of photons in output mode b. 
{. . .) denotes the quantum mechanical expectation value 
while the bars refer to an average over all realizations of 
disorder. The need for two averages is due to the fact 
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that quantum fluctuations and random multiple scatter- 
ing are both stochastic processes. C^^^/ is the quantum 
optical generalization of the classical intensity correlation 
function that has been intensely investigated in multiple 
scattering experiments [gI]. Note that the quantum cor- 
relations are induced despite the fact that the response 
of the multiple scattering medium is fully linear. Thus, 
a non-classical resource combined with linear optics is 
quite generally sufficient for many quantum optics ex- 
periments, and even enable quantum computing [l6| . 

The spatial quantum correlation function defined in 
Eq. ([T]) depends sensitively on the number and fluctua- 
tions of photons incident on the sample through direction 
a. Restricting to a theory for discrete spatial and tem- 
poral degrees of freedom, it is given by 
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{fla) is the average number of incident photons and the 
Fano factor Fa = ^n^al ("-a) gauges the variance in the 
photon number fluctuations specifled within the spectral 
bandwidth relevant for the measurement. C(^) and C(^) 
are classical spatial correlation functions Q that can be 
neglected in the diffusive regime of multiple scattering 
that is of concern here, but it would certainly be exciting 
to study their quantum optical equivalents as well. The 
size of the Fano factor determines the transition from the 
classical {Fa > 1) to the non-classical regime {Fa < 1). 
Using non-classical light opens the door to a genuine 

quantum regime where — 1 < C^haf < 0' corresponding 
to spatially anti-correlated photons. Such spatial quan- 
tum correlations are of infinite range in the sense that 
the magnitude is independent of the angular difference of 
the two output directions, which translates into a spatial 
separation in the far field. This is an example of the fun- 
damentally new phenomena that arise in quantum optical 
descriptions of multiple light scattering. So far, quantum 
optical variables have been measured in very few multi- 
ple scattering experiments that only have focused on the 
regime of classical photon fluctuations 12, 1^, while spa- 
tial quantum correlations were not demonstrated. 

The spatial quantum correlation of Eq. ^ can be ob- 
tained by recording the photon fluctuations of the total 
transmission An^. They are related through 



|R) + ^(n„)Vl,„ (3) 



which holds in the diffusive regime of multiple scattering 
where the spatial correlation function is independent of 
the output direction b. A similar equation applies to the 
case of total reflection measurements. Thus, the spatial 

quantum correlation function C^^^, is directly accessible 
by measuring in addition to the photon fluctuations also 
the average number of photons entering the sample {fia) 



and the ratio of the transport mean free path to the sam- 
ple thickness Note that the latter ratio is modified by 
so-called extrapolation lengths in order to account for the 
effects of the interface of the multiple scattering medium 
17| . In the experiment we extract the spatial correlation 
function in the frequency domain, which has dimensions 
of an inverse bandwidth. The complete continuous mode 
theory will be described elsewhere 
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FIG. 2: (color online), (a) Sketch of the experimental setup. 
Vacuum squeezing is generated in an optical parametric os- 
cillator with a PPKTP nonlinear crystal, overlapped with a 
displacement beam on a beam splitter (BS), and directed onto 
the sample (S). The multiple scattered light is imaged onto 
a photo-detector and its photon fluctuations are recorded us- 
ing an electronic spectrum analyzer (SA). A similar setup is 
used for reflection measurements illuminating the back side 
of the sample, (b) Measured photon fluctuations after mul- 
tiple scattering recorded in reflection geometry {L — 20 ^m) . 
The photon fluctuations of the non-classical light source (blue 
circles) are recorded relative to the photon fluctuations asso- 
ciated with the classical limit (gray triangles). Depending on 
the phase of a displacement beam, photon fluctuations below 
or above the classical limit are detected. The detection fre- 
quency was 3.93 MHz and the resolution bandwidth 300 kHz. 

As a non-classical resource we use squeezed light gen- 
erated in a second-order nonlinear process at a wave- 
length of A = 1064nm [l^. This versatile source features 
continuous tuning between classical (super-Poissonian) 
and non-classical (sub-Poissonian) photon fluctuations 
by varying the relative phase between the squeezed vac- 
uum state generated in the nonlinear process and a bright 
displacement beam. We are able to reduce photon fluc- 
tuations {Fa — 0.52 ± 0.02) below the classical limit 
{Fa = 1) and induce excess fluctuations above the clas- 
sical limit {Fa = 4.6 ± 0.2). The multiple scattering ex- 
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periment is conducted by focusing the non-classical light 
source onto either the front surface or the back surface 
of the sample, cf. Fig. [2Ia), to perform total transmis- 
sion and total reflection measurements, respectively. We 
record the total number of photons and the photon fluc- 
tuations transmitted through or reflected from the sam- 
ple. The multiple scattering samples consist of Ti02 that 
has been grounded into strongly scattering particles with 
a typical size of 200 nm. A range of samples with differ- 
ent thicknesses has been fabricated, and the total trans- 
mission of each sample was recorded with an integrating 
sphere. From these measurements the transport mean 
free path is found to be ^ = (0.9 ± 0.3) ^m and all sam- 
ples are in the diffusive regime where A/27r < £ ^ L. 

Figure Elb) shows the measured photon fluctuations, 
recorded in reflection, for squeezed light illumination of 
the multiple scattering medium. The classical limit was 
recorded by blocking the squeezed beam. The phase of 
the displacement beam is scanned in order to continu- 
ously tune the photon fluctuations of the incident light 
below and above the classical limit. We record photon 
fluctuations that are reduced 0.5 ±0.1 dB below the clas- 
sical limit after the photons have been multiple scattered, 
which corresponds to a Fano factor for the total reflection 
of Fjf = 0.90 ± 0.02. Correcting for the overall detection 
efficiency of (37 ± 1)%, we infer Fn = 0.72 ± 0.04 for 
light exiting the multiple scattering sample. This reduc- 
tion is the direct experimental proof that non-classical 
properties of light survive the complex stochastic process 
of multiple scattering even after ensemble averaging over 
realizations of disorder despite the common belief that 
quantum properties of light are fragile. We calculate the 
average number of scattering events in the multiple scat- 
tering process [20], which is given by ~ {L/iY . For the 
particular sample used for the measurements of Fig.[2l^b) 
we get N « 1000, including effects of the interface, which 
quantifies the complexity of the multiple scattering pro- 
cess. The observed reduction in the photon fluctuations 
shows that the information capacity associated with the 
multiple scattering channels can be enhanced beyond the 
classical limit, as predicted theoretically [l^ . 

We have carried out a detailed investigation of the 
transport of non-classical and classical photon fluctua- 
tions through the multiple scattering medium for a range 
of different sample thicknesses. Figure [3] displays the de- 
tected photon fluctuations after multiple scattering rel- 
ative to the fluctuations associated with the classical 
limit, which are plotted versus the sample thickness. Us- 
ing light sources with classical fluctuations, the multiple 
scattered light always displays excess photon fluctuations 
corresponding to the classical regime (Fig. [3ja)). Non- 
classical light allows entering the quantum regime where 
the photon fluctuations are reduced below the classical 
limit, see Fig. [3]Jb). Our experimental results can be 
compared to the predictions from the full quantum the- 
ory for multiple scattering of photons ll|, ll4| . Excellent 
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FIG. 3: (color online). Measured transmitted (red triangles) 
and reflected (blue circles) photon fluctuations after multiple 
scattering of light with classical (a) and non-classical (b) pho- 
ton fluctuations versus sample thickness. The data are pre- 
sented relative to the classical limit, which marks the bound- 
ary between the classical regime and the quantum regime. 
Every data point is obtained after ensemble averaging over 
six different sample positions. The black data points are com- 
parison to theory taking into account the collection efficiency 
of the detection setup, which was measured separately. The 
collection efficiency varies for each measurement point since it 
depends on sample thickness. The boundaries of the shaded 
areas correspond to the theoretical predictions for minimum 
and maximum efficiencies. 



agreement between experiment and theory is apparent 
from Fig. [3] both in the classical and quantum regime. It 
should be stressed that the comparison to theory requires 
no adjustable parameters, and only depends on the mea- 
sured photon fluctuations of the incident light source, the 
intensity transmission coefficient, and the total collection 
efficiency of the setup. 

The theory predicts that the multiple scattered pho- 
ton fluctuations depend on two terms varying as and 
L~^, respectively, cf. Eq. ([3]). The latter term is only ob- 
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lent agreement with theory. Controlling the power of the 
non-classical light source thus provides an efficient way of 
tuning the strength of the spatial quantum correlations. 
Furthermore, the spatial correlation function is predicted 
to be independent of the sample thickness, which holds in 
the diffusive regime of multiple scattering. This behavior 
is experimentally confirmed as well, see Fig. |4Ub). 

We presented the first experimental demonstration 
that photon fluctuations below the classical limit can be 
observed after multiple scattering using squeezed light 
as non-classical resource. Multiple scattering was found 
to induce infinite range spatial correlations that are of 
purely quantum origin. Our experimental results were in 
excellent agreement with the full quantum optics theory 
for multiple scattered light. We believe that our results 
will inspire further explorations of the multitude of new 
phenomena encoded in the quantum optical properties of 
multiple scattered light. 
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FIG. 4: (color online), (a) Measured spatial quantum correla- 
tion function C^^^, versus power of the incident light beam. 
For classical (black points) and non-classical (red triangles) 
photon fluctuations, positive and negative spatial correlations 
are observed, respectively. Every data point represents an av- 
erage over three different positions on the sample of thickness 
L = 6 fim. The curves are the theoretical predictions and 
the dashed line represents the uncorrelated case, (b) Spatial 
quantum correlation function versus sample thickness taken 
at an input power of P = 120 /^W. The spatial quantum 
correlation function is found to be independent of the sample 
thickness in agreement with theory (horizontal lines). 



servable for non vanishing spatial quantum correlations, 
i.e. the comparison between experiment and theory al- 
lows extracting the spatial quantum correlation function, 
cf. Fig. m We observe negative (positive) spatial corre- 
lations in the case where the transmitted photon fluctu- 
ations are in the quantum (classical) regime. According 
to Eq. ([2]) the strength of the spatial quantum correla- 
tion is expected to increase when reducing the number 
of photons of the incident light. This pronounced behav- 
ior is clearly demonstrated in Fig. [3]Ja) and is in excel- 
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